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Abstract. 

We suggest a direct algorithm for searching the Lax pairs for nonlinear integrable 
equations. It is effective for both continuous and discrete models. The first operator 
of the Lax pair corresponding to a given nonlinear equation is found immediately, 
coinciding with the linearization of the considered nonlinear equation. The second 
one is obtained as an invariant manifold to the linearized equation. A surprisingly 
simple relation between the second operator of the Lax pair and the recursion operator 
is discussed: the recursion operator can immediately be found from the Lax pair. 
Examples considered in the article are convincing evidence that the found Lax pairs 
differ from the classical ones. The examples also show that the suggested objects are 
true Lax pairs which allow the construction of infinite series of conservation laws and 
hierarchies of higher symmetries. In the case of the hyperbolic type partial differential 
equation our algorithm is slightly modified; in order to construct the Lax pairs from 
the invariant manifolds we use the cutting off conditions for the corresponding infinite 
Laplace sequence. The efficiency of the method is illustrated by application to some 
equations given in the Svinolupov-Sokolov classification list for which the Lax pairs 
and the recursion operators have not been found earlier. 
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In the present article we suggest an algorithm for constructing the Lax pairs for nonlinear 
integrable models. Our scheme is based on the symmetry approach. Let us explain the 
algorithm with the evolutionary type PDE of the form 

d^u 


Ut = f{x,t,U,Ui,U2, ...,Mfc), 
Recall that equation 

U-f tf 'll, Hi, 112, ..., Hm) 


dx^ 


( 1 . 1 ) 


( 1 . 2 ) 


is called a symmetry for the equation fll.ll) if function g satishes the following differential 
equation 


A- 


du dui 


A - ... - 


du, ‘ 


9 = 0 


(1,3) 


where Dt and are the operators of the total derivatives with respect to t and x 
correspondingly. Note that equation fll.Sp is overdetermined and for any hxed value of 
m one can effectively hnd all of its solutions of the form g = g[x, t, u, Ui, U2, (see 

[I])- 

An ordinary differential equation 

Um f^(^) ^5 Hi; ■■■I Hm—l) (^•^) 

dehnes an invariant manifold for fll.ll) if the following condition is satished 

AG-/^r/l(i.i),(i.4) = 0. (1.5) 

Here DtG is evaluated by means of the equation (II.Ih and all the x-derivatives of the 
order greater than m — 1 are expressed due to the equation fll.4p and its differential 
consequences. Equation (11.51) generates a PDE with unknown G admitting a large class 
of solutions. However it is a hard problem to hnd these solutions since the equation is 
not overdetermined. Some of the invariant manifolds to dm can be found from the 
stationary part of the symmetry fll.2p by taking g{x,t,u,Ui, = 0. 

We now concentrate on the linearization of the equation fll.ll) around its arbitrary 
solution u = u{x, t) : 


df df 


Hi -...- 


K 

duk 


Vk = 0 . 


( 1 . 6 ) 


Actually (11.61) dehnes a family of diherential equations, depending on u. The important 
fact that u ranges the whole set of the solutions to fll.ip is formalized in the following way. 
We assume that in addition to the natural independent dynamical variables h,Hi,H 2 , ... 
the variables u,Ui,U 2 ,... are also considered as independent ones. 

Let us hnd the invariant manifold of the form 


m—1 


Vm= 2_^ a{j)Vj, 
j=0 


V = Vo, 


(1.7) 
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to the equation fll.61) . This means that the condition 


' m—1 


D, 


d: 


^ ^ f 

I l(l-l)dl-6),(l-7) = 0 


( 1 . 8 ) 


Y1 

J=0 

holds identically for all values of u,Ui,U 2 ,.... Here we assume that Vj function 
a(j) depends on x,t and a hnite number of the dynamical variables u,Ui,U 2 , ■■■■ In 
(ll.Sp the variables Ut,Uxt,Uxxt, ■■■ are expressed by means of the equation ( 11 . 11 ) and 
the variables Vt, Vi^t ,■■■, are exresssed due to fll. 6 l) . Equation fll.81) splits 

down into a system of m partial differential equations with the unknown functions 
a(0), a(l),..., a{m — 1), coefficients of the decomposition (11.71) . Due to the presence of 
the additional dynamical variables u,Ui,U 2 ,... the system is overdetermined and hence 
can effectively be investigated. Actually, equation (11.71) dehnes a bundle of the manifolds 
depending on inhnite set of the variables u,Ui,U 2 , ■■■■ Suppose that such an invariant 
manifold is found. Then we can interpret a pair of the equations fll.611 . fll.711 as the Lax 
pair to the equation fll.ip . In the examples the order k of the equation fll.ll) and the 
order m of the equation fll.7p coincide. 

As an illustrative example we consider the well known KdV equation 

Ut = U'i + uui. (1.9) 

Linearized equation 

Vt = V'i + UVi + UiV (1-10) 

obtained by virtue of the rule (11.61) admits the third order invariant manifold dehned 
by the equation 


Vs = —V2 - + A 

Ml vs 


Vi + 


2 , “2 

-U + X] - Ml 

3 /Ml 


V. 


( 1 . 11 ) 


It is easily checked that equations fll.lOp and fll.lip are consistent if and only if 
the function m = u{x,t) satishes the equation fll.9p . Therefore these equations 
constitute the Lax pair to the KdV equation. It differs from the usual one found 
earlier in | 2 ]. Stress that there is no any second order invariant manifold of the form 
V 2 = a{u,ui,U 2 )vi + 6 (m,mi,M 2 )m for the equation (ll.lOp . But there is a hrst order one 
= f^v which however does not contain A and therefore does not generate any true 
Lax pair. 

Recursion operator is an important attribute of the integrability theory. It gives 
a compact description for the hierarchies of both symmetries and conservation laws. 
Various methods for studying the recursion operators can be found in literature (see, 
for instance 0-ini and the references therein). 

We observe that invariant manifold for the linearized equation is closely connected 
with the recursion operator for the original equation. Indeed examples in §6 show that 
equation dehning the invariant manifold, which provides the second operator of the Lax 
pair, can be rewritten as a formal eigenvalue problem of the form 


Rv = Am 































On a method for constructing the Lax pairs for nonlinear integrable eguations 


4 


for the recursion operator R. For instance, equation fll.lip is easily rewritten as (see 

§ 6 ) 

(Dl + 3 “ + 

where the operator at the l.h.s. is nothing else but the recursion operator for the KdV 
hierarchy. 

Therefore our scheme of constructing the Lax pairs provides an alternative tool 
for searching the recursion operator. On the other hand side when the recursion 
operator is known the invariant manifold and hence the Lax pair can be found by simple 
manipulations (see §6). At this point we have an intersection with the pioneering articles 
[1], where nonstandard Lax pairs are given for some of the KdV type equations in terms 
of the recursion operators. However these Lax pairs differ from those found within our 
scheme since they are nonlocal and do not contain any spectral parameter. 

There is a great variety of approaches for searching the Lax pairs from the Zakharov- 
Shabat dressing [121 US] and prolongation structures by Wahlquist and Estabrook [TT] to 
3D consistency approach developed in na-iiTi. We mention also approaches proposed 
in [iHiiiniii]. An advantage of our scheme is that it can be applied to any integrable 
model (at least in 1+1 -dimensional case), the hrst of the operators is easily found and 
the second is effectively computed. They allow Ending the conservation laws, higher 
symmetries and invariant surfaces for the corresponding nonlinear equation. The found 
Lax pairs are more complicated since they are based on the difierential operators of the 
orders greater than usual ones. This is their disadvantage. The question remains open 
whether the Lax pairs of this kind allow to find any new solution for the well studied 
models. 

In the case of the hyperbolic type integrable equations the algorithm should be 
slightly modified. Let us explain it with the example of the sine-Gordon equation 

Uxy = sinu. (1-12) 

Let us find the simplest but nontrivial, i.e. depending on a parameter, invariant manifold 
for the linearized equation 


Vxy = {cOSu)v. 


(1.13) 


It is the three dimensional surface in the space of the dynamical variables 
v,Vx,Vy,Vxx,Vyy, ■■■ defined by the following two linear equations, with the coefficients, 
depending on the field variable u{x,t) and its derivatives 


u 


Vyy — Uyicot U)Vy T X , Vx “ AU = 0 , 


smu 


Vxx — Ux{cot u)vx + A 


-1 


Un 


sin w 


Vy — X = 0. 


(1.14) 

(1.15) 


Here A is a complex parameter. Note that equations fll.l4p . fll.lSp are not independent. 
One of them is immediately found from the other by differentiation by means of the 
equations fll.l2p . fll.l3p . A triple of the equations fll.13p - fll.15p can be rewritten (see 
the end of §3) as a pair of the systems of ordinary differential equations providing the Lax 
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pair realized in 3 x 3 matrices. The method for deriving the Lax pair from the triple 
fll.13p - fll.15p is based on constructing the inhnite Laplace cascade for the linearized 
equation fll.l3j) and obtaining the finite reduction of the cascade. 

Let us give a brief comment on the structure of the article. In §§2,3 we discuss the 
well-known Laplace cascade for linear and nonlinear hyperbolic type equations. In §3 
also the problem of hnding hnite reductions of the inhnite Laplace sequence is studied. 
The Lax pair to the sine-Gordon equation is derived from the Laplace cascade. In §4 the 
dehnition of the invariant manifold for the hyperbolic type equations is recalled. The 
Lax pair is construted via invariant manifolds for hyperbolic equation fl4.12p found 
in ESI. In §5 the Lax pairs are constructed by evaluating invariant manifolds for 
the evolutionary type integrable equations. Beside the explanatory examples here we 
consider two equations fl5.19p and fl5.2np found in [21] as equations possessing inhnite 
hierarchies of conserved quantities. To the best of our knowledge the Lax pairs for the 
equations fl4.12p . fl5.19p and fl5.20p have never been found before. In §6 we illustrate 
applications of the newly found Lax pairs. The Lax pair obtained in the previous sections 
is used to construct conservation laws for a Volterra type chain. We also show that the 
second operators of our Lax pairs are closely connected with the recursion operators 
for the associated nonlinear equations. In Appendix we give all of the computational 
details appeared when we evaluated the invariant manifold for the linearization of the 
sine-Gordon equation. 


2. Laplace cascade for the linear hyperbolic type equations 


Let us recall the main steps of the Laplace cascade method (see [22], [23]). Gonsider a 
linear second order hyperbolic type PDE of the form 

v^y + a{x, y)v^ -h h{x, y)vy + c{x, y)v = 0. (2.1) 

It can easily be checked that functions 


^[ 0 ] = ttx + ah — c, /c[o] = by + ab — c (2.2) 

do not change under the linear transformation v —?• X{x,y)v with arbitrary smooth 
factors X{x,y) applied to equation fl2.ip . They are called the Laplace invariants for 

(El]). 

We rewrite equation fl2.ip as a system of two equations: 


d 




^[ 1 ] 


= hrniU. 


(2.3) 


When the invariant h[o] does not vanish then one can exclude v from fl2.3p and 
obtain a linear PDE for n[i] 

V[i]xy + + b[i]V[i]y + C[i]n[i] = 0 (2.4) 

where the coefficients are evaluated as follows 

d 


a[i] 


= a 


dy 


log(^[o]), b^i] = b, C[i] = a[i]6[i] b^i]y - h[o]. 


(2.5) 
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Thus we define a transformation of the equation fl2.ip into the equation fl2.4p . This 
transformation is called the Laplace ^/-transformation. Iterations of the transformation 
generate a sequence of the equations 


'^[i]xy T T CjjjUjj] 0 ( 2 . 6 ) 

for i > 1 where the coefficients are given by 

d 

a[i\ = ^ ^[i] = C[j] = a[j]6[i]-F6[i]j^-h[j_i].(2.7) 

Here we assume that a[o] = a, 6[o] = b, C[o] = c. Eigenfunctions U[j] are related by the 
equations 


d 

+ a[i] ] = U[i+i], 


9 , 




( 2 . 8 ) 


Due to the relation C[i] = ^a[i] + — h[j] system fl2.7p is rewritten as 


d 


a\i\ = a[j_i] — — log(h[j_i]), h[j] = h[j_i] -7 a[i\x — &[i] = b. (2.9) 

Reasonings above define the functions ajj], 6[j], h.[j] only for z > 1. However they can be 
prolonged for z < 0 by virtue of the same formulas rewritten as follows 

d 

h[i-i] = /Z[j] - a[i]a; - by, a[i_i] = — log(h[j_i]), 6[i_i] = &[,]. (2.10) 

Summarising the computations above we get a dynamical system of the form 
d 


dy 


log(^[i]) = ®[i] ~ ®[i+l ]5 a[i]* = ^[i] — ^[i- 1 ] + b[iT^=b, ( 2 . 11 ) 


which is reduced to the well-known Toda lattice 
d 

— log(h[i]) = p[i] - P[i+i], P[i]x = ^[z] - 

where p[i] = a[j] — 6 and b^ = by. 

Define two linear operators 


( 2 . 12 ) 


Li Dy Qj^i^ Dj. -|- 1 ] h^i_\^Dj^ , 


(2.13) 


where Dy are the operators of differentiation with respect to x, y correspondingly 
and Di is the shift operator acting as follows Dia\i\ = ajj+ij, Dj/if*] = h[j_|_i], etc. We 
summarize all the reasonings above as a statement. 

Proposition 1. The operators Lj, Mj commute for all z iff their coefficients satisfy 
the system fl^ . 

Corollary. Equations (12.81) constitute the Lax pair for the system (12.lip . 

The Laplace ^-transformation can be interpreted in a similar way. 
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3. Laplace cascade for the nonlinear hyperbolic type equations. 

Formal Lax pairs 


Let us explain how the Laplace cascade is adopted to nonlinear case [23] (see also [2l]h 
Consider a second order nonlinear hyperbolic type PDE 


Uxy F {xly I Uxi Urf). 


(3,1) 


Its linearization around a solution u{x,y) derived by substituting u = u{x,y,e) = 

on ^ y -j ) 

’ ’ |£=o into fl3.ip is an equation of the 


u{x^ y, 0) + ev{x, y) + ... with v{x, y) = 
form 


de 


Vxy + aVx + bvy + cn = 0 (3.2) 

where the coefficients a = —dF/dux, b = —dF/duy, c = —dF/du depend explicitly on 
the independent variables x, y and the dynamical variables u, Ux, Uy. Let us assign 
the Laplace sequence fl2.8l) - fl2.1ip to the linearized equation fl3.2p . However now instead 
of the operators d/dx, d/dy in (I2.8I) - (I2.11I) we use the operators D^i Dy of the total 
differentiation with respect to x, and y. Denote through Ui, Ui, i = 0,1,... the i-th 
order derivatives of the variable u with respect to x and y correspondingly 


DI.U = Ui, DyU = Ui- (3.3) 

Evidently we have explicit expressions for the operators Dx, Dy acting on the class of 
smooth functions of x,y and a hnite number of the dynamical variables Ui,Ui 


d 


d 


d 




i=0 


dui 


2=1 


dui 


(3.4) 


n _ ^ ^ a 

2=0 ^ 


CXD 


dy 


2=1 




(3.5) 


The Laplace invariants corresponding to the equation fl3.2p are evaluated as 

h[o] = Dx{a) + ab — c, fc[o] = Dy{b) + ab — c. (3.6) 

The linear system (12.8p in this case converts into 

{Dy + = n[j+i], {Dx + 6[j])n[i+i] = h[i\V[i\. (3.7) 

The coefficients aj*], 6[j], h[j] are evaluated due to the equations 

Dy{^og{h^if)') ci[j] Dx{a^if) /i[j] T Dy{bf 

b[t\ = b, a[o] = a. 

All of the mixed derivatives Ux;y,Uxxy, ■ ■ ■ are replaced by means of the equation fl3.1l) 
and its differential consequences. 

Inhnite-dimensional system fl3.7p dehnes a sequence of the linear operators 


Li Dy T ci[j] Di, Mi Dx T ij h^i_i'^Dj^ , 


( 3 . 9 ) 
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satisfying the commutativity conditions 

V* [Li,M,]=0. (3.10) 

Thus one can dehne a pair of commuting operators Lj, M* depending on an integer 
parameter i for an arbitrarily chosen equation fl3.ip . Roughly speaking the sequence 
of the commuting operators fl3.9p . fl3.10p recovers equation fl3.ip . Hence system fl3.7p 
dehnes a (formal) Lax pair for the arbitrary (generally non-integrable) equation (I3.ip . It 
is not very surprising since for non-integrable case system fl3.7p is of inhnite dimension. 
However as it is approved below by several examples for the integrable case the system 
is either hnite (Liouville type equations) or admits a hnite dimensional reduction (sine- 
Gordon type equations). 

Example 1. As an illustrative example of the non-integrable equation with the 
inhnite dimensional Lax pair consider the equation 

Uxy = U^. (3.11) 

For its linearization 

Vxy = 2uv (3-12) 

we have a[o] = &[o] = 0, h[o] = h[_i] = 2u, a[i] = — h[i] = u + One can hnd 

all of the coefficients ay], hy] due to the equations fl3.8p as functions of the dynamical 
variables and therefore dehne completely the system (13.7p . Now go back, suppose that 
the system evaluated above is consistent and show that its consistency dehnes uniquely 
equation (13.lip . Indeed the consistency of (13.7p implies Dx^ay]) = hy] — h]o] equivalent 
to - u which gives (13.111) . 

Example 2. As an example with the hnite system (13.7p we take the Liouville 
equation 

Uxy = e“ (3.13) 

for which h]o] = h]_i] = e“ and hy] = /i[_ 2 ] = 0. For i > 1 and i < —2 the Laplace 
invariants hy] are not dehned. The coefficient ay] is dehned only for the following three 
values of i : aqj = —Uy, a[o] = 0, a[_i] = Uy. Evidently b]Q] = 0. System (13.7p for the 
equation (I3.13p contains only seven equations 

{Dy - Uy)vy] = V] 2 ], DxV] 2 ] = hy]Vy] 

DyV]o] = n[i], DxVy] = h]o]V]o] 

{Dy + Uy)v]-i] = n[o], DxV]o] = h]_i]V]-i] 

DxV]-l] = h]_2]V]_2] 

There is a freedom in choosing of vyy Put them equal to zero. Then the 

obtained system gives the Lax pair for (I3.13P 

T, = AT, Tj, = RT, 

where T = (n[i],n[o],and 
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/ 

0 


0 



( 

Uy 

0 

0 

\ 


0 

0 



,B = 


1 

0 

0 


V 

0 

0 

0 

) 



0 

1 

— Uy 

) 


A = 


Remark that there are some degenerate cases where the commutativity condition of 
the operators fl3.9p dehnes not exactly the initial equation fl3.ip but some other equation 
connected with (13.ip by a Miura type transformation. Illustrate it with the following 
example. 

Example 3. Consider the equation 




(3.14) 


The coefficients of its linearization 


-'xy 


= e’'(v + Vr) 


(3.15) 

assigned 

to fl3.14p implies the equation r^y = e^ + r^^, connected with fl3.14p by a very simple 
Miura type transformation Ux + u = r. 


depend on r = u + Ux- Thus the commutativity condition of the operators 


(3.16) 


3.1. A more symmetrical form of the Laplace seguence 

Let us change the dependent variables in the system fl3.7l) to make formulas more 
symmetrical. Introduce new dependent variables i G (—cxo, cx)) in such a way 
that tC[j] = V[i] for z > 0 and W[j] = h[_i]h[_ 2 ]...h[j]z;[j] for i < —1. 

Then the set of equations fl3.7p . (13.Sp is changed to the form below, where z > 0 : 

{Dy + a[i\)wii] = W[i+i], {Dx + 6[o])'W^[j+i] = h[i\wii], 

Dy{\og{hii])) = a[i] — a[j+i], Dx{a[i]) = h[j] — /z[j_i] + T)y(6[o]), 

{Dx + {Dy + a[o])'u;[_i_i] = h[-i-i]W[-i], 

= 6[_i] - Dx{\og{h[_i_i])), 6[o] = ^[0] 

3.2. Sine-Gordon eguation 

More than two decades ago an important property of the Laplace invariants of the 
Liouville type integrable equations has been observed PI, pg. It was proved that 
the hyperbolic equation (13.11) is an integrable equation of the Liouville type if and 
only if the set of its Laplace invariants is terminated on both sides. Mention also 
recent results obtained in [26] . The problem of describing the properties of the Laplace 
invariants characterizing the sine-Gordon type integrable PDE is discussed in [27]. Our 
investigation convinces that a connection between the sine-Gordon type equations and 
the Laplace cascade is clearly formulated in terms of the cascade eigenfunctions. Let us 
explain our observation with an example. 

We consider the sine-Gordon equation 


Uxy = smzz. 


(3.17) 
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It can be shown that the Laplace invariants h[j] for the linearized eqnation 

Vxy = {cosu)v (3.18) 

do not vanish identically for any integer i. Tims system (13 .Th provides an inhnite- 


dimensional Lax pair for fl3.17p . Below we show that for this case fl3.7p admits a 
hnite-dimensional redaction. Bring fl3.7p to the symmetric form fl3.16p . 

Proposition 2. The system (I3.16p corresponding to the sine-Gordon eqnation 
(13.17p with tCfo] = a[o] = 0, = 0, C[o] = —cosu is consistent with the following 

cntting off bonndary conditions 

W[2] = a(-l)M;[_i] + a(0)M;[o] + a(l)tf^[i], 

W[_2] = /?(-l)w[-i] + /3(0)w;[o] + /?(1 )m;[i] 


(3.19) 


where 


COS u smu ^ 


«(-l) = «(0) = A, a(l) = 


(3.20) 


A is a complex parameter. 

Sketch of proof. Look for the fnnctions (13.201) providing the consistency of the 
following overdetermined system of eqnations obtained from (13.161) by imposing (I3.19P 


'^[i]y = («(1) - + «(0)w[o] + a(-l)tf^[-i], 

W[o]y = wp], 

W[-l]y = h[-l]W[ 0 ], 


(3.21) 


and 


(3.22) 


W[i]x = h[o]-u;[o], 

W[0]x = 

W[_i]x = + /?(0)m;[o] + (/?(-!) - 6[_i])m;[_i]. 

The compatibility conditions {w[i]x)y = {w[i]y)x generate a system of nonlinear eqnations 
for the fnnctions a(j),/3(j) searched: 

D,(/3(-l))+/3(l)«(-l) = h[_2], 

Il,(/3(0)) + /3(l)«(0)+/?(-l)h[_i] = 0, 

Dy{P{l)) + /3(l)(a(l) - ap]) +/3(0) = 0, 

Dx{ot{—l)) + «(—1)(/3(—1 ) — &[-!]) + q;( 0) = 0, 

Dx{oi{0)) + a;(—1)/3(0) + a(l)h[o] = 0, 

L>a;(a(l)) + a(-l)/?(l) = h[i]. 


(3.23) 


Here all the given coefficients apj, &[-i], h^- 2 ], ^[o]) ^[i] sxe linear fnnctions of the 

derivatives Ux, Uy\ 


h[o\ = h[_i] = cosM, h[i] = h[_ 2 ] = — 
0(1]= My tan M, 6 [_i] = Ma-tanw, 


+ 


COR^ 'll, 5 


COS U 


(3.24) 
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therefore we can assume that C({j),/3{j) also linearly depend on the hrst derivatives of 
u 

a{j) = u, Uy) = p{ 3 , u)uy + g(j, u), ^ 5 ) 

/3{j) =/3{j,u,u^) = r{j,u)u^ +s{j,u), j = l,0,-1. 

Substitute expressions (13.24^ . (13.251) into equations (13.231) and then compare the 
coefficients before the independent combinations of the dynamical variables u^Uy, Ux,Uy. 
As a result one gets twenty four equations for twelve functions p{j,u), q{j,u), r{j,u), 
sUiU), j = 1,0,—1, depending on u only. We write down explicitly only a part of 
the equations since the others are obtained from these by applying the replacement 
p{j) O r(-j), q{j) O s{-j). 


!S:^ + r(l)p(-l) = jy^, + s(l)p(-l) = 0. 

= 0, r(-l) sinu + s(l)g(-l) = 

«) + r(l)p(0) = 0 , « + s(l)p(0) = 0, 

r(—1 ) cosM + r(l)g(0) = 0, r(0) sinn + s(—1 ) cosu + s(l)g(0) = 0, 
^ - r(l) tann + r(l)p(l) = 0, ^ - s(l) tanu + s(l)p(l) = 0, 

r(0) + r(l)g(l) = 0, r(l) sinu + s(0) + s(l)q'(l) = 0. 


(3.26) 


By solving the overdetermined system of equations (I3.26p we hnd explicit expressions 
fl3.20p for the coefficients of the constraint fl3.19p . Now the systems fl3.2ip . fl3.22p can 
be rewritten as follows 


= 5 ^ 

where 4/ = (tcq], tcjo], W[_i])^ and 


/ 


A = 


0 

0 


\ A sir 


COSM 

0 

1 

X 


0 

1 

Ut cot u 




/ 


,5 = 


/ 


V 


(3.27) 


Uy cot U 
1 
0 


A 

0 

COSM 


sin u 

0 

0 




/ 


. (3.28) 


It is easily checked that (13.271) . (I3.28p dehnes the Lax pair for the sine-Gordon equation 
fl3.17p . We failed to reduce it to the well-known usual one found in 


4. Invariant manifolds of the hyperbolic type PDE 

Let us recall the dehnition of the invariant manifold of the hyperbolic type equation 
(13.ip . Consider an equation of the form 

G{x, y, Uk, Uk-l, ...U, Ml, M2, • • • Urn) = 0. (4.1) 

Note that G depends on x,y and a set of the dynamical variables u,ui,ui ,..., where 
Uj = 1^, Uj = 1^. Take the differential consequences of fl4.ll) 

Gi{x, y, Uk+l, ...u, Ml, M 2 , ... Urn) = 0 , 

G 2 {x, y, Uk,... M, Ml, M2, . . . M^+l) = 0, 


(4.2) 

(4.3) 
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where Gi,G 2 are evaluated by applying the operators D^^Dy : Gi = 0 ^ 0,02 = DyG 
and subsequent replacement of the mixed derivatives by means of the equation fl3.1l) 
and its differential consequences. Equation (14.Ih dehnes an invariant manifold for (13.11) 
if the following equation is satished 


E*a;E*i;G'|(3.1),(4.1)-(4.3) — 0. 


(4.4) 


Example 4. Show that equation 


iixx 2 tan xi — 0 


(4.5) 


dehnes an invariant manifold for the sine-Gordon equation fl3.17l) . Here 
G = Mxx + lu^tanu, Gi = D^G = + \ul and G 2 = DyG = ^ + 

easily verihed that 


2(cos u) 


It is 


DxDyG 


Dx 



ulUy \ 
2(cosm)^/ 


0 mod((4.5),G = 0,Gi = 0,^2 = 0). 


Therefore equation fl4.4p holds and thus fl4.5p dehnes an invariant manifold for fl3.17p . 


4 . 1 . From the Laplace cascade to invariant manifolds 

Show that reduced system fl3.16p . fl3.19p is closely connected with the invariant 
manifolds of the linearized equation fl3.18p . Indeed, equations fl3.16p imply that 
W[2] = {Dy + a[i]){Dy + a[0])'u;[0], W[i] = {Dy + a[o])'u;[o], W[_i] = {Dx + 6[o])'u;[o], 
W[_ 2 ] = {Dx+b^-i]){Dx+b^)w[o]. Therefore since a[o] = &[o] = 0, the boundary conditions 
(I3.19P turn into the equations 

{Dy + aii])DyWio] = a{-l)DxWio] + a(0)w[o] + a{l)DyWio], 

{Dx + bi_i])DxW[o] = (3{-l)DxWio] + /3(0)w[o] + (3{l)DyWio]. 

Simplify the equations obtained by using explicit expressions fl3.20p . fl3.24p and hnd 
equations 

u 

Ly'^lO] ■ Uy cot uDy T A—: Dx A)'U7[o] I (3.17),ui[Q]j.y = (cosn)ui[Q] 0, (4-6) 

SlU u 
XI 

Lx'^[0] • X tix cot -|- A —: Dy A )'W^[0] I (3.17),'UI[0]3:y —(‘^GS ^ (^*^) 

siri xjj 

which dehne the invariant manifold for fl3.18l) discussed in Introduction (see fll.Mp . 
fll.lSp above). 

This observation leads to an alternative algorithm to look for the Lax pair. Instead 
of the cutting oh boundary conditions to the lattice fl3.16p one searches an invariant 
manifold for the linearized equation (I3.18p . 

By construction we have 

Mw[o] := {DxDy - cosu)w[o]. 


(4.8) 
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Commutators of the operators M satisfy the following relations 

[-hx) -hj/] ‘2Axy(^X Ly \Lx)i 

[M, Lx] = BxxLy — BxyLx + {Axx — A ^Axy)M, (4-9) 

[M, Ly] = ByyLx “ BxyLy + {Ayy “ XAxy)M 

where A = logcot|, B = logsinu, Ax = Dx(A), Ay = Dy(A), Axx = B>^(A) and so 

on. Consequently any element of the Lie ring generated by the operators Lx, Ly, M is 
represented as a linear combination of the same three operators. 

Linear equations (14.bh . (14.7p dehne a manifold parametrized by ^[o], W[Q\y and 
the dynamical variables u, Ui, Ui, .... By applying Dx to the equations fl4.6l) . (14.71) and 
then simplifying due to the equations (13.171) . (I3.18p one gets another parametrization of 
the manifold 

'^[0]xxx "^[Olxx “ 1 “ ^x)'^[0]x ^ ^[ 0 ]: 

Ux Ux 

Sin u sin u 

'^[0]y ^ “1“ ^(cOS'u)'U^[o]a; ^[0]? 

^tJjx ^^x 

where the parameters tC[o], tC[o] 3 ;, are taken as independent ones. It is shown below 

that this parametrization is closely connected with the Lax pair for the potential KdV 
equation being a symmetry of the sine-Gordon equation. 

4-2. Evaluation of the invariant manifolds and the Lax pair for the eguation 
Uxy = f{u)^,yi + ul, f" = 7/ 

In this section we construct a Lax pair to the equation 

Uxy = f{u)^l + ul, f" = 7/ (4.12) 

found in [20]. It is known that S-integrable equation of the form (I4.12p by an appropriate 
point transformation can be reduced either to the case f{u) = u or f{u) = sinu (see 
12®. By analogy with the sine-Gordon equation considered in the previous section we 
look for the invariant manifold of the form 

Vyy + aVy + Ijvx A cv = Q (4-13) 

for the linearized equation 

Uxy = f\u) \/l + ulv + -^^^Vx. (4.14) 

Apply the operator Dx to (I4.13P and rewrite the result as 

Vxx = {2vUxf{u)f'{u) + Vxf{u) + Dx{h)vx + cVx + Dx{c)v + Dx{a)vy) 

{uyUxVx + {Vy + av){l + ul))f'{u) + {pUyv{l + ul) + aUxVx)f{u) 

b./TT^x ^ 

Now apply Dy to (14.151) . simplify the result due to the equations above and get an 
equation of the form 

Vyy + dVy + bVx + CV = 0, 


(4.10) 

(4.11) 


(4.16) 
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with the coefficients a, b, c depending on a hnite number of the dynamical variables. 
According to the dehnition of the invariant manifold equations fl4.13p and fl4.16p should 
coincide. This fact implies a system of three equations on the sought functions a, b, c 

(2f{u)uxbf{u) + D^{c)b - Dy{h)D^{a) + DyD^{a)b - abD^{a)^ a/TT^ 

+ {2puybul - D^{a)u^b + 2puyb)f{u) - {ul + l)Dy{h)f'{u) = 0, (4.17) 

{i.-Dy{h)ul + ah- Dy{h))f{uf + UyhiZ + 2ul) f {u) f {u)^ + 

-(b^D^ia) - Dy{c)b - DyD,{h)h + Dy{h)c + Dy{h)D,{h)){l + ulf/^ 
Dy{jX^Ux^ Dy{h^(X'llx^{\ ~\~ 

+Ux{l + ul){b‘^Ux + Uyyb + aUyb - Dy{b)uy)f{u) = 0, (4.18) 

(2pf{uyuyUxb + Ux{ab - 2Dy{b))f{u)f{u) + 3f{u)‘^Ua:Uyb 

-Dy{b)D^{c) - chD^{a) + DyD^{c)h^ + ul + h{ul + 3)f'{u)f{uf 

+ {{-Dy{h)a + Dy{a)h + D^{h)h + pulh){l + ul) + f\u) 

+ {p{buyy + aUyb + - Dy{h)uy){l + - D ^{c)uj)) f{u) = 0. (4.19) 

Assuming that the searched functions depend only on u,Ux, Uy i.e. a = a{u, u^i Uy), 
b = h{u,Ux,Uy) and c = c{u,Ux,Uy) substitute these functions into (I4.17p . (I4.18p and 
(I4.19P and eliminate mixed derivatives of u using (I4.12p from the resulting equations. 
Thus we obtain three equations of the following form 

Otiiu^ Uxi Uy)UxxUyy T (di{Ui Uxi Uy)Uxx T '^y)'^yy T diiu, Ux-, Uy) 0, 

i = 1, 2, 3. These relations are satished only if the following conditions: 

ai(u, Ux, Uy) = 0, fti{u, Ux, Uy) = 0, 7i(M, Ux, Uy) = 0, 6i{u, Ux, Uy) = 0 (4.20) 

hold identically for all values u, Ux and Uy, i = 1,2, 3. Here 

Oil {boii^uy OjxiJjuy) \/l T ul, 

02 = {bbu^uy - 

O 3 {bOu^Uy Ou^buy)-\/l + ul, 

/^i T bUyOi^n^ bn^UyOn^ aboy^^) ■\l 1 T Ux 

+ f{u){l + ul)(bay^y^ - by^OuJ, 

1^2 (It ^x)^^u'Uybux bUybyy^ T 6 T 

+ f{u) {{bby^y^ - blj{l + ul)^ + bby^Ux{l + ul) + b^) , 

1^3 fbUyCyy^ Cbtty^ byUyCy^)ljl T 

+ (1 + uDibCy^y^ - by^CyJfiu) T 6 (6^, (1 T M^) T bUx)f{u), 

71 Ux(^ Oiubuy T ^®tillj;)\/l T 'U'x 

+ (1 + ul){{hay^y^ - bu^au^)f{u) - f{u)by^), 

72 = (1 + ul) 


{-buybuUx + bbuuyUx - f{ufbuy - cbuy + 6c„Ja/1 + ul 
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+ {{bhuyUy aU^Ky + UMuy)f{u) 

+ u^{b-buyUy)f{u) , 

73 = Ux{-Cubuy - 2buyf{u)f{u) + bcuuy)V^ + “x 

+ (1 + Ul){{pb - buyCuy + bCuyUy - PU yb U y) f (u) 

+ {-abuy + bauy)f'{u)), 


5i = 


(1 + ul){-bu,auy + bau^uy)f{uy + {-u%, + 2u^b - bujf{u)f{u) 

'^xiPu'^y^U abOy^ bCy bUyayy^ \/ 1 “ 1 “ uy, 

+ b{ ^y^uuy + Cyy + 2pUy — aOyy + ^UUx^X )il + nl)f{u) 
bay - {b 

Ux + byHyOy ){1 + U 


U 


So = 


+ (1 + ul){bay^u^ - byUy + bay^Uy)f'{u), 

ibh 




huXy){l + <Yf\u) 


+ (1 + ul){bay^Ux - aby^U^ + bbyyUa, - byUy)f{u) 

+ abpiu) + Uy{ 2 ulb + 3 b - ulby^ - Uxbyjf{u)f{u) 

-(1 + ul){-bUybyyU^ + + b^UyU^ + CbyUy - bCyUy) \/1 + 

-|- (1 -|" ^x) PbCy^ “1“ bby b ttyy bbyyyUy bybyyUy T" bbyy^Ux 

- cby^ - by^byUx “ by^f{u)){l + U^) + U^UylbOy - Oby + pbUy) 


S. = 


+ {{bby^Uy + bby^Ux){l ^ u^) + + UxaUyb - UxbyUy)f{u) 

(1 + ul){{bCy^yy - pUyby^ - by^Cyy) + 2 pUyUxb) f (u) 

+ {{boy^ + bbyy - abyj{l + ul) - 2uxbyUy + Uxab)f'{u)f{u) 

~\~ 3f [vfjbUxUy Uxi^cbay byCytly bCyytiy'j \/1 “ 1 “ 

+ (36 + bul - 2by^ul - 2by^Ux)f{u)f{u) 

+ ((1 + ul){pUxb^ + bUyCyyy - C 6 a„j, - byUyCyy - by^CyUx “ PU^y 

+bcyy^Ux + apUyb) + bcy) f {u) 


- (1 + ul){abyUy - bCy^Ux - bbyUx - bOyUy - bCyyUy - pbUy)f{u). 

Thus the problem is reduced to a system of equations fl 4 . 20 p . 

We look for the functions a, b and c depending on the variable Uy linearly: 

a a\{u^ Ux)Uy T 02 (n, Uxf) , 

6 b\{xi^ Ux^Uy 62 (n, Ux ^5 
C = Ci{u,Ux)Uy + C2{u,Ux). 

Then equation 0:2 = 0 is essentially simplified 62(6i)„^ — 6i(62)„^ = 0 . Assume that 
62 = 0, then 

b = bi{u,Ux)uy. ( 4 . 21 ) 
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- c„, = 0 . 


OjUx 0 ) 

Assume that = Cu^ = 0 then 

a = ai{u)uy + a 2 {u), c = Ci{u)uy + C 2 {u). (4.22) 

Substituting the functions fl4.2ip . fl4.22|) into 71 = 0 (see fl4.20l above)) we obtain 
the following equation 

bi{u,u^) {f{u) + ai{u)f{u)){l + M^) + a/1 + uluo^a'^iu) ) = 0. 

Since functions /, oi and 02 depend only on u, we obtain 
f'{u) + ai{u)f{u) = 0 , ^(m) = 0 . 

Consequently 

/ ^ /'(“) / N 

ai[u) = — 77 ^, a 2 (M) = 03 , 

fw 

where 03 is an arbitrary constant. Then from the equality 73 = 0 we get 
ci(u) = C 2 {u) = -f{u) + C 3 , 

fw 

where C 3 is an arbitrary constant. Analyzing the equation 5i = 0 we dehne that 
b — ^3 

where 63 7 ^ 0 is an arbitrary constant. From the equation ^2 = 0 we obtain that 03 = 0 
and C 3 = — 63 . It is easily checked that equalities (3i = 0, i = 1, 2, 3, 72 = 0, ^3 = 0 are 
automatically satished. 

Thus summarizing the reasonings above we can claim that equations 


■'yy 


f{u) 






f{u)^Jl + ul 


{f{u) + X)v = 0, 


, f (^) I Uxx 
'^xx I ^ “r 




Ux\/ul + 1 


fiu) ' ul + lj— ' Xfiu) V A 

dehne an invariant manifold for the linearized equation (14.141) . Here A is the spectral 
parameter. 

Derive the Lax pair for the equation (I4.12p from the invariant manifold (I4.23p . 
fl4.24p . To this end evaluate the Laplace sequence of the form fl3.16p for the linearized 
equation (I4.14p . In what follows we will need in explicit expressions for several hrst 
coefficients of the system fl3.16l) 

Uxfiu) 


ul + ^ 


(4.23) 


u = 0 (4.24) 


®[o] = ®[-i] = 




&[o] = b[i] = 6[_i] = 0, 


7 -eU ^ ! 2 i 1 V, '^xXf{u) 

^[0] = / (m)v“x + m = -——-T + 




(u: 


«[i] = 


+ l)f{u){'yuy^/ul + l- f{u)) + UxUxxP{u) - f'{u)UyUxx\/ul + l 

^/ul + l{f'{u){ul + 1) - f{u)Uxx) 


(4.25) 
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The constraints fl4.23p . fl4.24p generate the cntting off bonndary conditions of the form 

np] = a(l)n[i] + a(0)n[o] + a(-l)n[_i], 

V[-2] = /3(l)w[i] + /3(0)n[o] + /3(-l)r;[_i] 

imposed on the inhnite system fl3.16p . Find the coefficients a{i) and fi{i) in the relation 
(I4.26p . Evidently eqnations (I3.16p imply 

{Dy + a[i\){Dy + a[o])n[o] = np], {D^ + + 6[o])^^[o] = ^[-2]- (4.27) 

Set V := np] and simplify (I4.27p by virtne of (I4.26p and the relations n[i] = {Dy + a[o])n 
and n[_i] = {D^ + 6 [o])n. As a resnlt we obtain 

Vyy + (ap] + np] — Q!(l))ny — a(—l)^^; + (apj^ + np]Ci[i] — ap]Q!(l) — a(0) — a(—l)6p])n = ( 
Vxx + (&[-!] - ld{-l))vx - ld{l)vy + {b[o]x + b[o]b[_i] - ap]/3(l) - /?(0) - (3{-l)b[o])v = 0, 
The last eqnations shonld coincide with fl4.23p . fl4.24p . Comparison of the corresponding 
coefficients allows one to derive explicit formnlas for the songht fnnctions a{i) and (3{i) 

a(l) = ap] + ap] + 
a(-l) = 


\Uy 


/(n)\/ 7 |+T’ 

Q!(0) = f{uY + X +ay- ah - 


+ 


/W f{u)y/ul + l ’ 


c\/ul + l 


m = 

I3{-1) = 6[_i] + 


DO) 


/3(0) = b[o]x - a[o]( 


Uxf'ju) _ 

f(u) y ui+ i ’ 


)-( 


Uxf'ju) UxU, 


Yb[0] + 


ul+l 


\f(u) ! '' f(u) ' ul+l mVJ] I A 

Now the inhnite system of eqnations 03.161) is rednced to a pair of the third order systems 
of ordinary differential eqnations 

V[o]y = '^[ 1 ] ~ «[o]i^[o], 

Di]y = («(1) - «[i])^[i] + «(0)^^[o] + 

V[-i]y = fcpjnp] - ap]n[_i], 

V[o]x = '^’[- 1 ] - ^[i]'y[i], 

■f^plx = ^[o]'i^[o] - ^[o]'y[i], 

n[_i]^ = /3(l)^^p] + /3(0)np] + (/?(-!) - 

which can be specihed as follows 

/ Uxfju) ^ 

%/'W 


^ ^[0] ^ 


^[1] 


V ^1-4 ) 

y 

^ ^[0] ^ 


V[i] 

= 


aAI +1 

A 


0 




V 


_ Uxfju) 

fO) a/7|+T 
f{u)^yul + l 0 


Uxfju) 

Y “^+1 / 



(4.29) 


/ 0 

fin) 


V ^[-1] / 


\AI+i 


Uxxfju) 

oi+iY 


0 

0 


1 

0 


Ux^ul + l Uxf'ju) I UxUxx 




f{u) ul + 1 


\ ( 
) 


^[ 0 ] ^ 
^[ 1 ] 

V ^[-11 / 


(4.30) 


’ (4.28) 


Systems 04.291) . 04.3Up dehne the Lax pair for the eqnation 04.121) . 
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5. Searching the Lax pairs for the evolutionary type integrable equations 

Let us consider the evolutionary type integrable equations, for which the Laplace cascade 
is not dehned. Here we use the scheme set out in the Introduction. 

5.1. Korteweg-de Vries eguation. 

As an illustrative example we consider the Korteweg-de Vries equation 

Uf "^xxx ~h VjVjx. (5.1) 

Its linearization evidently has the form 

Vxxx "h UVx "1“ UxV. (5.2) 

Direct computations show that equation 05.21) does not admit any invariant manifold of 
the form Vxx = ci{u,Ux,Uxx)vx + b{u,Ux,Uxx)v ht for arbitrary solution u{x,t) of (15.Ih . 
Let us look for the invariant manifold of order three 

Vxxx = aVxx + bvx + cn, (5.3) 

where the coefficients a, b, c depend on a hnite number of the dynamical variables u, 
Ml, U 2 , ■■■ ■ According to the dehnition the following condition 

{Vxxx}t ~ {yt)xxx (^•^) 

should be valid. 

Replacing in 05.4p Vxxx and Vt due to 05.3p and 05.2p respectively and then comparing 
the coefficients before the independent variables Vxx-, Vx and v we obtain the following 
equations 

3aDx{b) + Quxx + Dl{a) + 3Dx{a)b + UxO -|- uDx{a) + 3Dx{a)^ 

+3aD‘l.{a) -|- 3ofDx{a) + 3Dx{c) — Dt{a) + 3Dl(b) = 0, (5.5) 

^xib) + 3D‘^(c) + 3bDx{b) -\- 3Dx{a)Dx{b) -|- uDx{b) + ^Uxxx + 3abDx{a) 

+3D‘l.{a)b — 3aUxx — Dt{b) + 3cDx{a) + 2uxb = 0, (5.6) 

3acL)x(^a'j -|- Uxxxx nOxi^cj -I- L)^(^cj -I- 3L)nuxxx A 3L)x(Jf)c 

-Dt{c) + 3uxC + 3Dx{a)Dx{c) - buxx = 0. (5.7) 

It is reasonable to assume that a = a(u,Ux,Uxx), b = b{u,Ux,Uxx) and c = c{u,Ux,Uxx)- 
We substitute these expressions into 05.51) . 05.61) and 05.71) and then exclude all the mixed 
derivatives of u due to the equation 05. ip . As a result we obtain three equations of the 
following form 

Q(j('U, Mx, Uxx) l^xxx^^xxxx jdiiu., ^xx')^xxx V(^) ^xi ^xx)'^xxx 

diiu., Uxy Uxx'j'^xxx ^xi '^xx) 0, i 1; 2, 3. 

Since Uxxxx, u^xxxi '^‘xxx) '^xxx are independent variables then these equations split down 
into hfteen equations as follows 

Ux^ Uxxi Uxxx) 0) (diiUi UxT Uxx) 0, Uxj Uxx) 0, 

dii^UjUxjUxx') 0, Ci( m, Mj,,'Ua;^) 0 


(5.8) 
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hold for all values of u, and u^x, i = 1,2,3. Thus the searched coefficients a, 6, 
c satisfy a highly overdetermined system of differential equations fIS.Sp . Specify and 
analyse the system. It can be verified that the three equations /3i = 0, i = 1,2,3 
immediately imply = 0, = 0 and = 0 therefore 

a Ux)Uxx T U2(u, Ux)Vjxx 0>‘^{Vj^ 1 

h = hi{u,Ux)ul^ + b2{u,Ux)Uxx + b3{u,Ux), 

C Cl (m, Ux)u^^ C2 (Uf Ux)Uxx T C3 (u, Ux) ■ 

Equations a, = 0, * = 1,2, 3 are of the form 

^UxxUxx'^XXX + 0,0,Uxx T ^UUxx'^X T T O'UxUxx'^XX = 0, (5-9) 

buxxUxx'^XXX T Cuxx T buUxx'^X T buxUxx'^XX T ^0,Uxx ~ 0, (5.10) 

‘iCuxxUxx'^XXX + ^CuUxx'^X T ‘^CuxUxx'^XX T ^('^Uxx +1 = 0. (5.11) 

Then since functions a, b and c depend only on the variables u, Ux and Uxx the coefficients 
at Uxxx vanish, i.e. we get 


a Ui('U, Vjx)Uxx + U2(u, Ux) •; 
b bi (m, Ux)Uxx + ^2 (^, Ux) 1 
C Cl (u, Ux)Uxx + C2 (u, Ux) ■ 

Substituting a, b and c into the equations (15.Oh . (15.101) and (I5.1ip we obtain 

{al + {ai)ux)uxx + + ^ 3 ,( 01 )^ + &! = 0 , 

((^i)^^: + (^ibi)uxx + Cl + Ux{bi)u + a,ib2 = 0, 

3((ci),i^ + ciiCi)**^^ + ‘d>Ux{ci)u + 3aiC2 + 1 = 0. 


Since functions Oj, bi and c,, * = 1,2 depend only on u and Ux the coefficients at Uxx 
vanish and we get the following system of equations 

+ {(^i)ux = 0, 

{bi)ux + = 0, 

(c.)„. + a.Ci = 0 , 

0x02 + Ux{ai)u + 5i = 0, 

Cl + Ux{bi)u + ciib2 = 0 = 0, 

Sux{ci)u + 3aiC2 + 1 = 0. 


Concentrate on the last system. It is easy to check that ai 7^ 0. Assume that &i 7^ 0 
and Cl 7^ 0. From system fl5.12l) we find 


ai = 


u. 


-\- a 


u] 


0,2 — 


Ur 


(Cl)n + bi 


bi = 

b2 = 


bi{u) 


Ux + 04 (u) ’ 

Cl + Ux(bl) u 


Cl = 


C2 = 


C4(m) 


Ux + 04 ( 1 *) ’ 

1 + ^Ux{Ci),^ 


(5.13) 


Ui Ui 3cii 

Now equations 7^ = 0, * = 1,2,3 are satisfied automatically. From equation ^i = 0 we 
obtain that 64 = 65 and 04 = a^u + a^, where 65, 05 and qq are arbitrary constants. 
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From equation ^2 = 0 we get 04 = |m + C 5 , where C 5 is an arbitrary constant. Equation 
^3 = 0 implies 05 = 65 . Then from equation ei = 0 we obtain 05 = Oe = 65 = 0. It is 
easy to check that equations 62 = 0 and £3 = 0 are identically satished. 

Thus equation fl5.3p is of the form 


Uxx ( 2 

= - Vxx - 0“ + -^ 

Ux. V o 


2 , \ Ux^x 

-u + x) — 


Un 


(5.14) 


Here A = C 5 is an arbitrary parameter. 

Proposition 3. Pair of eguations h5.X\) . \5.1f\) defines the Lax pair for the KdV 
eguation. 


5.2. Potential and modified KdV eguations 
Concentrate on the potential KdV equation 


^xxx T 2 ^^ 


Its linearization 


Vt = Vxxx + where w = Ux 


(5.15) 


(6.16) 


does not admit any second order invariant manifold of the necessary form Vxx = 
a{u, Ux, Uxx, ■■■)vx + h{u, Ux, Uxx, However it admits a third order invariant manifold 
given by 

Wx / o , \ ,UJx 

w 


Wx / 9 , \ ^Wx 

= —Vxx - {w + A) + A— V. 
w w 


(5.17) 


Here A is an arbitrary parameter. The consistency condition of the equations (Ih.lbh . 
(15.17p is equivalent to the mKdV equation 

1 


Wt Wxxx "^ 2 ^ 


(6.18) 


connected with fl5.15p by a very simple substitution w = Ux- In other words fl5.16p . 
fl5.17p dehne the Lax pair for the equation fl5.18p as well. 

Now let us return to the sine-Gordon equation fl3.17p . Recall that equation fl5.15p is 
a symmetry of the sine-Gordon equation. It is easily seen that equation fl5.17p coincides 
with (H.lUp up to the notations. 


5.3. Lax pairs for the KdV type eguations from Svinolupov-Sokolov list 


Gonsider the following two third order differential equations 

3 

^yyy 7,'^y 


Ut = Uyyy “ —Uy ~ '^f VPUy, 

^UxUlx 7 3 

Ur = Uxxx -;-:TT- Ux 

2{l+ ul) 2 


(5.19) 

(5.20) 
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possessing infinite hierarchies of conservation laws [21] • As it is established in [20] these 
equations are symmetries of the equation fl4.12l) . We have proved above in §4 that 
equations 

f'{u) Am,, 


^yy ft ( \ ^y'^y 


fiu 


f{u)y/l + ul 


Vx - if (u) + A)m = 0, 




f'{u) u^x 


+ 




Urr,\/Ul + 1 


Vy - 


X 


(5.21) 


M = 0 (5.22) 


f{u) ul + lj Xf{u) 

dehne an invariant manifold for the linearized equation (I4.14p . It is reasonable to expect 
that invariant manifolds for the linearizations 

Vt = Vyyy - ^(qW^ + f{u))Vy - 3/(m)/'( 

and 

3 ( (l- ul)nl 

2 V (1 + ulT- 


3UxUxx / V-^ '^xJ'^xx , 2 1 

I^T Vxxx T ^ y "Vxx 7> 1 7i i ~9T9 a X'^x I 


1 + M 


(5.23) 

(6.24) 


of the symmetries (I5.19p and (I5.2np are closely connected with the same manifold. Indeed 
by applying the operators Dy and Dx to the equations fl5.2ip and respectively to (15.221) 
one can deduce the following two third order ordinary differential equations 


u 


-'yyy 


^Vyy - {^U^y + nu) + X)Vy 


u 


{P{u) + \)Uyy 


M,. 


3/(m)/'(m)% P = 0, 


(5.25) 


(1 + 3ul)uxx 

Vx.x.x. / . X Vx 


{l + ul)Ux 


2„,2 


, 1 9 UxUxxx ^uiu 

A"^ +jul + 


u 


, -^xx 

Vx H-A M = 0. 


(5.26) 


1 + ul {l + ulYj Ux 
It is easily checked by a direct computation that equations fl5.25p and (15.261) dehne 
invariant manifolds for (I5.23p and (I5.24p correspondingly. Conclude the reasonings with 
the following statement. 

Proposition 4. 1) Linear eguations / 1 5. 23\) . 15. ^51) define the Lax pair for the 
eguation 

2) linear eguations \5.2Jif), h5.2dt) define the Lax pair for the eguation H5.20\) . 


5.4- Volterra type integrable chains 

In this section we discuss the semi-discrete equations of the form -^Un = 
f{un+i,Un,Un-i) with the sought function u = Un{t), depending on the discrete n and 
continuous t. The direct method for constructing the Lax pairs through linearization 
can be applied to the discrete models as well. As illustrative examples we consider the 
modihed Volterra chain 
dpn 


dt 


= -P„(Pn+l -Pn-l) 


(5.27) 
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and the equation 
dUn 


(5.28) 


dt Uji—i 

found in [29]. These two equations are related to each other by a very simple Miura 
type transformation 


dur. 


= Pn 


Pn = 


dt Un—X 

Note that the coefficients of the linearization 


dvn 

dt 


= -Pl{vn+1 -Vn-l) 


(5.29) 


(5.30) 


of the equation (I5.28p depend on the variable pn- This explains why we study these two 
equations together. Look for the invariant manifold of the third order: 

Vn+2 = aVn+1 + bVn + CVn-1 (5.31) 

to the equation (15.301) with the coefficients a, b, c depending on a hnite set of the 
dynamical variables Pn,Pn±i, ■■■ ■ Actually we suppose that fl5.3ip dehnes an invariant 
manifold for any choice of the solution p = Pn(t) to the equation (15.271) . 

The coefficients a, b, c are found from the equation 
d 


dt 


{aVn+l + bVn + CVn-l) = DI [-pl{Vn+l - Vn-l)) ■ 


(5.32) 


Studying the equation fl5.32p we assume that the variables {Pfc}^_oo 5 Vn+i^ Vn-i are 
independent dynamical variables. Omitting the simple but tediously long computations 
we give only the answer: 

A , . A 


a = 


Pr 


+ 


Pn+l Pn+1 PnPn+l 

Therefore the invariant manifold searched is of the form: 
Pn A 


b = 1 


c = 


Pn 

Pn+1 


Vn+2 — 


Pn+1 


+ 


pI+1 


1 - 


Vn+1 

A 


Vn + 


Pn 


1 • 


PnPn-\-l / Pn-\-l 

Proposition 5. The consistency condition of the eguations H5.30\) and (E 
coincides with the eguation |(5.^7D . 


(5.33) 


6. Construction of the recursion operators and conservation laws via newly 
found Lax pairs 


It was observed that the Lax pairs for the integrable equations found above essentially 
differ from their classical counterparts. In this section we discuss some useful properties 
of the newly found Lax pairs. We show, for instance, that they provide a very convenient 
tool for searching the recursion operators and conservation laws for integrable models. 
As illustrative examples we take the KdV and pKdV equations, the Volterra type chain 
fl5.28p etc. We show that the equation of the invariant manifold to the linearized equation 
is easily transformed into the recursion operator. 
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Let us start with the KdV equation (11.91) . We can rewrite equation (ll.ll|l of the invariant 
manifold in the following form 


Uxx 2m ^ 2 uUx 

Dl-—Dl + —Dx + Ux-^ 
Ur 3 iUr 


V = XUxDx- 


-V. 


tier 


( 6 . 1 ) 


We now multiply 


from the left by the operator UxD^^— and obtain a formal 


eigenvalue problem of the form 
Rv = Xv 
for the operator 


R = (dI- ^Dl + 

Ux \ Ux 


2 uux 

3Ur 


( 6 . 2 ) 


(6.3) 


An amazing fact is that R coincides with the recursion operator for the KdV equation. 
Indeed, we have 


Rv = UxD~^ 


Vxx + UxD., ^ ^ + V - - 


— 1 / ^XXX^X 


'\Jxxx'^-^X %Jjxx'^'7 


ui 


+ - 


2 UVr 


3 Ur 


+ v - - 


2 UUrrV 


3 ui 


‘^xJ-'x 


2 UVr 


2 UUrrV 


3 Ux 3 
It can be simplihed due to the relation 


UVr UUxxV / UV 

u 


Ur 


ut 


and reduced to the form 


Rv= [ D^ + -u + -UxD^ 


-1 


V. 


In a similar way we derive from fl5.17p the recursion operator to the potential KdV 
equation (15.151) . Indeed rewrite (I5.17P as follows 

1 


u 


Ri - —DI + uiDx 1 V = —XuxD. 


Ur 


2 

x-^x 


Ur 


(6.4) 


Equation (16.4p implies 
Rv = —Xv 

where 


R = UxD:^— { Dl - —Dl + ulDx 




11:1 


Simplify the expression for Rv: 


Rv = UxD~^ 


— 1 / ^XXX^X 


'tJxxx'^^X ^Ujxx'^ 'Y 


Ut 


UxVx 1 Uxx UxD VjxDxV 


^ X ^x^x^ 


Apparently operator R = + UxD^ ^UxDx coincides with the recursion operator for 

dSISI). 
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Invariant manifolds fl5.25p and 05.261) for the linearized equations 05.23P and 05.241) 
allow constracting of the recursion operators 

Piu) 


R = 


tin, 


+ UyDy {jUyy-ff) 


and 


R = Dl-‘^?^ + u,Dp 


1 + ul 


'^XXX ^ ^x'^xx 


l + uO(l + nl)^ 
for the Svinolupov-Sokolov equations 05.19P and 05.20p respectively. 


x I L)x 


6.2. Recursion operators via the Lax pair for a Volterra type chain 


We proceed with an example of the discrete model 05.28p . Let us write equation 05.33p 
defining the invariant manifold for 05.30p as follows 



Pn 

Pn+l 


Dr. 


- 1 - 


Pn 

Pn+l 



Vn 


Pn+l 


{Dn-l) " 


Pr. 


(6.5) 


where Dn is the shift operator acting due to the rule Dna{n) = a{n + 1). We multiply 
06.5p from the left by the factor {Dn — l)~^Pn+i and then get 


Rv = Xv (6.6) 

where R = Pn {Dn — l)~^Pn+i (^D"^ + ■:^^Dn — 1 — ■^pf^Dn^'^ . After some elementary 
transformations the operator R is reduced to the form 

R = pI {Dn + Dp) + 2pnPn-l + 2pn (A>„ - 1)"^ {Pn-1 “ Pn+l) (6.7) 


Operator R given by 06.7p dehnes the recursion operator for the chain 05.28p . For 
instance, by applying the operator R to --we obtain the r.h.s. of the equation 

_ "fin + l— 'tLn — 1 

constructed in [30] 


—u - 1 ( 1 ^ 1 \ 

dr {Un+l — Un—l) \^n+2 Un Un—2 ) 


( 6 . 8 ) 


as a fifth-point symmetry of the equation 05.281) . Note that the whole hierarchy of the 
symmetries for the chain 05.28P is described in [3T] . 


6.3. Conservation laws via the Lax pair for a Volterra type chain 


Let us consider an equation of the form 
dun 1 

dt Un+l Un—1 


(6.9) 


Due to the relations 05.29P between u and p a pair of the equations O5.30|) . 05.33P define 
a Lax pair to the chain 06. 9 p as well. 

Rewrite the scalar Lax pair 05.301) . 05.33P in the matrix form 


Vn+l fVn^ 


dpn 

dt 


QUn: 


( 6 . 10 ) 
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where 


( 


f = 


g = 


p7i I A 
Pn + l P^+l 


1 - 


Pn 


1 

VO 10/ 

/ PnPn+l - A -PnPn+1 \ 


PnPn+1 Pn+l 

0 0 
1 




-Pn 


y Pn—lPn 


Pn 
0 

Pn-1 

Pn 


pI 


( 6 . 11 ) 


( 6 . 12 ) 


A -pn-lPn + A / 


To construct the conservation laws, we apply the method of the formal diagonalization 
suggested in [32l |33] and developed in |3l] . 

The first equation in (I6.10p has singular point X = oo {f has a pole at A = cx)). It 
can be checked that the potential / is represented as follows 

/ = 


where 


( ^ 


a = 


/5 = 


. 0 0 \ 

pt+i P"+i 

A-i 


V 

/1 
0 
VO 


0 

Pn+1 

Pn 

1 

0 


Pn + l g 
Pn 

1 


(6.13) 


1 / 


pnPn+1 

A—PnPn + 1 

Pn 

A PnPn-\-l 




/ 


(6.14) 


A — PnPn-\-l 

are analytic and non-degenerate around A = oo, Z is a diagonal matrix of the form 
/ A 0 0 \ 


Z = 


V 


0 1 0 
0 0 A-1 


(6.15) 


The change of variables tp = fdy is reducing the first system in fl6.1Up to the special form 

pn+l = PZpn, (6.16) 


where 


P = = 


( 


\ 


1 _ Pn+Pn+2 y-1 

Pn+1 P"+l 

A-1 

0 


Pn+2 

Pn 

Pn+1 

Pn 


1 - 
1 - 


PnPn+1 

A—Pn+iPn + 2 
PnPn+1 \ 


A—Pn + iPn+2 / 


Pn + 1 


A 


A—Pn+iPn + 2 


Pn+lPn+2 ^ 
A—Pn + iPn+2 
Pn + 1 

A—Pn+iPn + 2 
Pn+lA 
A—Pn + iPn+2 / 


The function P(A) together with P ^(A) are analytic around A = oo: 

OO 

P(A) = ^P^lX-\ 

where 


(6.17) 


*>0 


( ^ 


p(0) — 


Pn + 2 
Pn 
Pn+l 

Pn 
^2 


- 

Pn+1 

0 0 

.,2 


0 \ 


p(l) — 


( - P"+Ap^+P"+ 2) Pn+lPn+2 \ 

Pn + 1 

1 


V 0 / 


-Pn+1 


-Pn+1 


V 


0 


Pn+lPn+2 Pn+lPn+2 / 
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p(2) _ 

/ 

0 

P^+lP^+2 

pI+ipI+2 

\ 

, ^ 

/ 

0 

P^+lP^+2 

pI+iPI+2 

\ 



0 

-p^+lPn+2 

-pi+lPn+2 



0 

-p^+lP^+2 

-pi+lPl+2 


etc. 


V 

0 

Pn+lPn+2 

Pn+lPl+2 



V 

0 

Pl+lPn+2 

pI+ipI+2 

) 



The leading principal minors of -P(A) 

det P{X = oo) = det PiX = oo) = —^—, 

1 pI+1 2 

det P{X = oo) = det P{X = oo) = 

3 Pn 

do not vanish if the variable pn satisfy the inequality Pn 7 ^ 0 for all n. According to 
the Proposition 1 in [33] the hrst system in fl6.10p can be diagonalized, i.e. there exist 
formal series 

T = T® + + • • •, (6.18) 

h = + ■ • ■ (6.19) 

such that the formal change of the variables = Tip converts the system (I6.16p to the 
system of the diagonal form 

Pn+l = hZpn- ( 6 . 20 ) 

Thus we see that the formal change of variables y = Rp = reduces the hrst 

system in investigated Lax pair fl6.10p to the form fl6.20p . By construction R = f3~^T is 
a formal series of the form 


R = + R^^'>X-^ + ^ 

It follows from (I6.19I) - (I6.21I) that diagonalizable system (16.201) admits an asymptotic 
representation of the solution to the direct scattering problem 

yn{X) = R{n, A)e^"=~o ^°sHs,x)^n ^g_ 22 ) 


with the “amplitude” A = R{n, A) and “phase” 0 = nlogZ + A). 


s=no 


Turn back to the problem of diagonalization. By solving the following equation 

DniT)h = P(A)T, f = ZTZ-^ (6.23) 

we hnd the formal series T and h 


T = 



P„+lP'n+2 

Pn 

0 


0 PnPn-l{Pn - Pn-2 

Pl+lPl + 2iPn + l+Pn + 3) 


0 ^ 
0 

0 / 


A 


-1 


Pn 

0 


0 




4 9 

PnPn-l 


T 


( 2 ) 


32 


PnPn+1 

0 


A-2 + • • • 


/ 


h = 


'P — 

Pn + 1 


0 

0 


0 

Pn + 1 

Pn 

0 


0 \ 


pI+1 J 


( 


+ 


Pn+l{Pn+Pn+2) 

P^l 

0 

0 


0 

P'i+liPn + 2-Pn) 
Pn 

0 


0 

0 




A 


-1 


P^+l(Pn +Pn+2) / 
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/ 


V 


PnPn+2 

0 

0 


k 


( 2 ) 

22 


0 

0 


A-" + 


0 pt+li^PnPn+2 + pI + Pl+ 2 ) 


where 


T^2 = PnPn-l{pl- 2 Pl-l + Pl-2Pn-lPn-3 + pI-iPI + ‘^Pl-lPnPn-2 “ Pn+2pl+lPn), 

^( 2 ) pl+x{-pl+2Pn+3 - Pl+2Pn+l + pIPu-I + PnPn+lPn+ 2 ) 

^22 ~ 


Pr, 


According to the general scheme the second system of the Lax pair 06.101) is diagonalized 
by the same linear change of the variables y = Rip, where 


R = (3-^T = 

^ 0 
+ -pI 

\ Pn-lPl 


/ 

1 

Pn^l 

0 \ 



Pn 



0 

1 

0 


V 

0 

Pn—lPn 

Pn / 





Pn+lPn+2 


Pn 

PnPn-\-l 


PnPn-\-l 

-pI 


\ 


n-lP^Pn+1 - 

pI-iA ■ 

- Pn-2Pl-lP 

di’ 

p(2) 

-^12 

r'S \ 


p(2) 

-^21 

TjO) 

2^22 

p(2) 

^23 

A-2 + ---, 

p(2) 

\ R3I 

p(2) 

-^32 

flg’ ) 



A 


-1 


(6.24) 


(6.25) 


(6.26) 


?( 2 ) 

Li 

?( 2 ) _ 
h2 

?( 2 ) 


= PnPl+lPn+2, 

Pn+lPn+2{Pn+2Pn+3 ~h Pn+lPn+2 ~h PnPn+l) 


-^12 ~ 


Pn 


R\3 =PnPl+l{Pn + Pn+2), 

R?1 = -Pn(2Pn-l +Pn+l), 


(2) V V V V 

R22 = PnPn+lPn+2 + 3pn-lPnPn+l + PnPn+1^ 

R?3 = -^PlPn+1, 

R^il = 2pl-lPn + Pn-lPnPn+1 + Pn-2Pl-lPt 

R 32 = -Pn-lPn{^Pn-lPlPn+l + pIpI+1 + Pn-2Pn-lPnPn+l + pI-2 PI-1 
+Pn-3Pl-2Pn-l + pI-iPI + ‘^Pn-2pl-lPn) ■ 

This change of the variables reduces the second system in fl6.10p to the form ^ = Sy^ 
with 


S = -R-^Rt + R-^gR = 


-1. 


/ 1 0 0 \ 
0 0 0 
0 0-1 


A 



/ 

2p„p„+i 

0 

0 ^ 


d? 

0 

0 

\ 

+ 


0 

Pn{Pn—l Pn+l) 

0 

+ 

0 

c(2) 

‘-’22 

0 



V 

0 

0 

^PnPn-\-l j 


^ 0 

0 

CO 

CO to 

/ 


A-1 + - 
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*^11 PnPn+liPn+lPn+2 Pn—lPn)i 
*5*22 PnPn+liPn—lPn Pn+lPn+ 2 } ^ 

“S'!? = -PnPn+l{Pn+lPn+2 + Pn-lPn)- 
According the paper [33] the equation 
Dt\nh = {Dr,-l)S 

generates the inhnite series of conservation laws for the equations fl5.27l) . We write down 
in an explicit form three conservation laws from the infinite sequence obtained by the 
diagonalization procedure 
1 


A In 


Pn+l 


^')PnPn+li 


Dti-Pn+liPn + Pn+ 2 )) = {Dn “ ^)PnPn+liPn-lPn + Pn+lPn+ 2 ), 


Dt \^-2pnPn+lPn+2 “ -PnPn+1 “ ^Pn+lPn+2 

= {Dn - l)plpl^-i^{2pn-lPn+2 + Pn+lPn+2 + Pn-lPn)- 
Being rewritten in terms of these relations give the conservation laws for the equation 

dnn): 

1 

Df ln('U^_|_2 Un) {Dn 1) 


(^n+1 l)(^n+2 ^n) 


D, 


^n—1 ^n+S 


(^n+l l)(^n+2 ^n)(^n+3 ^n+l) 


= {Dn - 1) 


{Un Un—2'){Un+l ^n—l) (^n+2 ^n) 

1 


D, 


{'^n+l ^n—l)(^n+2 ^n) (^n+3 ^n+l)/ 

2 


(^n+l l)(^n+2 ^n) (^n+3 ^n+l) 

1 

Un—l) (Wn+2 ^n) 

1 

2('nn+2 ^n)^(^n+3 ^n+l)^ 

= {Dn - 1) 


{Un Un—2'){Un+l '^n—l) {'^n+2 '^n) (^n+3 ^n+l) 

1 


(^n+1 l)^(^n+2 ^n)^(^n+3 ^n+l) 

1 


+ 


"Uri,— 2 ) (^n+1 ^n—l)^(^n+2 ^n) / 

These conservation laws coincide with those found earlier (see, for instance, [18], [29]). 
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Conclusions 

There is a large set of classification methods allowing classes of integrable nonlinear 
PDEs and their discrete analognes to be described. For studying the analytical 
properties of these equations one needs the Lax pairs. Therefore the problem of creating 
convenient algorithms for constructing the Lax pairs is relevant. In the present article 
such a method is suggested. For the evolution type integrable equation the Lax pair 
consists of the linearized equation and the equation of its invariant manifold. In the 
case of the hyperbolic equations to obtain the Lax pair we use the Laplace cascade in 
addition to the invariant manifold. The method is applied to equations (I4.12p . (I5.19p 
and fl5.20p known to be integrable for which the Lax pairs have not been constructed 
before. 

An interesting observation is connected with the Laplace cascade of the sine- 
Gordon equation. It is proved that in this case the cascade admits a hnite-dimensional 
reduction which generates the Lax pair to the sine-Gordon model. We conjecture that 
the Laplace cascade corresponding to any hyperbolic type integrable equation admits a 
hnite-dimensional reduction. 

We considered examples showing that our method leads to true Lax pairs having 
useful applications. For the Lax pair of the Volterra type chain we found an asymptotic 
eigenfunction which allowed an inhnite set of conservation laws to be constructed. It 
is also shown that these Lax pairs allow the recursion operators to be constructed 
describing inhnite hierarchies of the higher symmetries and invariant manifolds for the 
given nonlinear equation. 
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7. Appendix. How to look for the invariant manifolds for the linearization 
of a nonlinear hyperbolic type equation 

It is well known that the linear hyperbolic type equation v^y = p{x,y)vx + q{x,y)vy -f 
r{x,y)v might admit linear invariant manifolds (see [35]). Here we consider equation 
fl3.2p obtained by linearizing an essentially nonlinear equation fl3.ip . We request that 
at least one of the coefficients a, b, c in (13. 2 p depends on at least one of the dynamical 
variables u, Ui, Ui. We look for a linear invariant manifold of the form 

{p{j)Dl+p{j-l)Dl-^+...+p{l)Dy+p{0)+q{l)D^+...+q{k)D^^)v = 0.(7.1) 

It is assumed that all the coefficients p{i) and q{m) can depend on x, y and on a hnite 
set of the dynamical variables u,Ui,Ui,U 2 ,U 2 , More precisely we study a family 
of equations (13.2p . (17.ip depending on a functional parameter u{x,y) such that when 
u = u{x,y) ranges a set of all solutions to the equation (13.ip then (17.ip ranges a set 


















On a method for constructing the Lax pairs for nonlinear integrable eguations 


30 


of invariant manifolds for the corresponding equation fl3.2p . In this case equation fl7.ip 
generates an overdetermined system of differential equations for dehning the searched 
coefficients p{i) and q{m). 

As an illustrative example we consider equation fl3.18p obtained by linearizing the 
sine-Gordon equation. We search an invariant manifold for (I3.18p in the following form 

Vyy + avy + Ijv x + cv = 0. (7.2) 

Here a, b and c are functions depending on a hnite number of the dynamical variables 
u,ui,ui,U 2 ,U 2 , ■■■■ Apply the operator Dx to (17.2p and rewrite the obtained result as 
follows 

Vxx = {sm{u)uyV - cos{u){vy + av) - Dx{a)vy - Dx{b)vx - Dx{c)v - cvx) ■ (7.3) 

Now apply the operator Dy to this equation and simplify the result by means of 
the equation (I3.18p . We arrive at the equation 

Vyy + dvy A bvx A cu = 0 , 

which should coincide according to the dehnition above with equation fl7.2l) . i.e. 
a = a, b = b, c = c. 

These conditions give rise to the following equalities 

DyDx{a)b — Dy{b) cos(m) — 2 sm{u)uyb A Dx{c)b — Dx{a){ab A Dy{b)) = 0, (7.4) 

DyDx{h)b - Dy{h)Dx{h) A Dy{c)h - Dy{b)c - h'^Dx{a) = 0, (7.5) 

DyDx{c)b — Dy{b)Dx{c) — cbDx{a) A cos(m) {—Uyb A Dy{a)b A Dx{b)b — Dy{b)a) 

A sin(-u) {—Uyyb — b'^Ux — auyb A Dy{b)uy) = 0. (7.6) 

Assume that a = a{u,Ux,Uy), b = b{u,Ux,Uy) and c = c{u,Ux,Uy) and substitute 
these functions into (17.41) . (17.5p and (17.6p . Eliminating the mixed derivatives of u due 
to equation fl3.17p we obtain three equations of the following form 

Uxj Vjy^Uxx^yy A ’^x'i '^y)’^xx ^xj ^y^'^yy 4“ di{u^ "^xj ^^) b? 

i = 1 , 2 , 3. Since the functions /3i, 7 * and 6i, i = 1,2,3 depend only on u, Ux and Uy, 
we should have the coefficients at UxxUyy, Uxx, Uyy and remaining terms equal to zero, 

i.e. 

Oli{U,Ux,Uy^ 0, (di{U,Ux,Uy^ 0, Vjx, Uy^ 0, Si{U,Ux,Uy^ 0 ^7.7^ 


for all u, Ux and Uy, i = 1,2, 3. Here 

ai = bou^uy buyOu^c, (7-8) 

^2 = bb^^uy - buybn,,, (7.9) 

0^3 bCu^Uy CuJOuy, (7.10) 

fdi = bcu,, - abau,, - bnUyOu,, - sin(M)a„^ 

A sin(M) A bauu,,Uy, (7.11) 

f ^2 = bbuu^Uy A bbu^u^ sin(M) - b^Uybu^ - b'^Ou,, - bl^ sm{u), (7.12) 
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/^s = sin(M) + hcuu^Uy + bbu, cos(m) - 

- bu, sin(M)c,,^ - buUyCu^, (7.13) 

71 = - cos{u)buy + bsm{u)auyuy + bu^Ouny - buyauU^ - buyOuy sm{u), 

72 = bUxbuuy - buybuUx + bsm{u)buyUy - bl^ sm{u) + bcuy - cb^y, 

73 = - buyCuUx; - bsm{u) - b^yCuy sm{u) + bouy cos(m) + sm{u)uybuy 

+ bUxCuuy - a cos{u)buy + bsm{u)cuyUy, 

= {by,^auy - ba^^uy) cos^(m) - cos(m) sin(M) 

+ {bOuyUy + bOu^Ux; - buUy) cos(m) 

+ {bou - 2uyb + bouuyUy - buUyauy) sin(M) 

“t“ ( Ojbauy by^^OyUx "1“ bcy^y^ sin.('u) 

“1“ buxauu'^y by^UyauUx nbay^Ux ba^^uy bu^nuy L~ bc^Ux-, 

S 2 = {bu^buy - bbu^uy) cos^(m) + b{bu,Ux + buyUy) cos(m) 

+ {-buUybuy - bu^buUx + buxbuu^) sin(M) 

+ {bbuuyUy - b'^auy - cbu, + bbu + bcu,) sin(M) 

“t“ bby^^y^y bCyUy b^UyUx 0) yU y 6 Oyllx “t“ bUxbyyUy by^by^^ 

bs = {by^Cuy - bCy^yy - UybyJ cos^(m) + {boy^ - aby^ + bbyy) cos(m) sin(M) 
+ ibbyUx - abyUy + boyUy + bcy^Ux + bcyyUy - u^b) cos(m) 

“t“ (b'^x^uux aUyb “1“ bcyyyUy by^CyUx) sin('?/) 

+ {-cbayy - byUyCyy - b^Ux + bcy + u^Oy) sin(M) 

“1“ bUxCyyUy byUyCyUx “t“ '^ybux by^Cyy “|“ bCy^yy . 

Hence the problem of searching the equation (17.2p is reduced to the system of equations 

(rTTD . 

We now look for the functions a, b and c depending linearly on the variable Uy: 

a = ai{u,Ux)uy + a 2 {u,Ux), (7.14a) 

b = bi{u,Ux)uy + b2{u,Ux), (7.146) 

C = Ci(u,Ux)Uy +C 2 iu,Ux). (7.14c) 


Then equation 02 = 0 (see fl7.9p i takes the form 62 (^ 1 )^^ — bi{b 2 )yx = 0. This equation 
is satisfied only if at least one of the following conditions holds 


61 = 0 , 62 

Let us consider the case 62 
the form 


( 6 l)n^ ( 62)1 


0 ’ h h ■ 

bi 02 

0. Then function b defined by the formula fl7.146p takes 


b = bi{u,Ux)uy. (7.15) 

Furthermore equations ai = 0 and eta = 0 become 

'^yO‘UxUx OjUx 0) '^y^UxUx 0- 
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These equations hold if at least one of the following four conditions is satisfied: 

= 0 , Cu, = 0 ; 

= 0 , Cu^ 7 ^ 0 , 

= 0 , 7 ^ 0 , 

4)0^3.7^ 0, 


^UxUx 

1 


Uy 

^UxUx 

1 _ 

0>Ux 

Uy 

OlUxUx 

^UxUx 




It can be proved that cases 2)-4) lead to contradiction. Concentrate on the case 1) 
which implies 

a = ai{u)uy + a 2 {u), c = ci{u)uy + C 2 {u). (7-16) 

By substituting the functions (17.151) . (I7.16p into equalities (3i = 0, i = 1,2,3 (see 
fl7.lip . fl7.121) and fl7.13p i one can verify that equation /3i = 0 is satished. The equations 
/92 = 0 and /ds = 0 lead to 

bi{bi)uu,, - {bi)u{bi)u,, = 0 and (cosm)m^ 6 i( 6 i)„^ = 0 

correspondingly. Thus bi{u,Ux) = Fi{u) for some function Fi{u). The equality 71 = 0 
becomes 

Fi(n)(cosM + UxO^iu) + ai(n) sinu) = 0 . 

Clearly, the functions ai and 02 are defined by the formulas 
aiiu) = —cotu, 02 (u) = Cl, 

where Ci is an arbitrary constant. According to that we rewrite the equality 73 = 0 as 
F 1 {u){uxc! 2 {u) + Ci(m) sinu + cosuCi) = 0 . 

Since u and Ux are regarded as independent variables the last equation leads to 
Ci(m) = —Cicotu, C 2 (m) = C 2 , 

where C 2 is an arbitrary constant. Let us turn to the equation (5i = 0 which gives: 

UyUx{Fi{u) cosu + F[{u) sinu) =0. 

Integration of the equation leads to 

^3 


Fi = - 
smu 

Now the equation ^2 = 0 takes the form 

C 3 Uy{CiUySmu + (C 3 + C 2 ) cosu) = 0 

We get C 2 + C 3 = 0, Cl = 0. It is easy to check that equalities ^3 = 0 and 72 = 0 are 
identically satished. Thus we have 


a = —UyCotu, b = X 


u 


y 


smu 


c = —A. 


Therefore the invariant manifold (17.2p is of the form 

u 

sinu 

and coincides with fl4.6p while fl7.3p gives fl4.7p . 


Vyy — Uy COt UVy + A . ^ “ AU = 0 
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